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Recently it has been demonstrated that a careful treatment of both longitudinal and transverse 
matrix elements in electron energy loss spectra can explain the mystery of relativistic effects on the 
magic angle. Here we show that there is an additional correction of order (Za) 2 where Z is the 
atomic number and a the fine structure constant, which is not necessarily small for heavy elements. 
Moreover, we suggest that macroscopic electrodynamic effects can give further corrections which 
can break the sample-independence of the magic angle. 



I. INTRODUCTION 

The title of this article is in reference to a recent work 
by Jouffrey et alX with the title "The Magic Angle: A 
Solved Mystery." The magic angle in electron energy 
loss spectroscopy (EELS) is a special value of the mi- 
croscope collection-angle a c at which the measured spec- 
trum "magically" becomes independent of the angle be- 
tween the incoming beam and the sample "c-axis." The 
mystery, in the context of 200 keV electron microscopy, 
is that standard semi-relativistic quantum theory yields 
a ratio of the magic angle 9m to "characteristic angle" 
Be of more than twice the observed 2 - value. Unfortu- 
nately, time^ and again,— the theoretical justification 
of the factor of two turned out to be an errant factor 
of two elsewhere in the calculation. A key contribution 
of Jouffrey et al. was the observation that relativistic 
"transverse" effects, when properly included in the the- 
ory, naturally give a factor of two correction to the non- 
relativistic magic angle. Here we show that there are yet 
additional corrections to the theory which can even break 
the sample independence of the magic angle. 

As in Ref. we consider here the problem of a 
relativistic probe electron scattering off of a macro- 
scopic condensed matter sample. Similar problems have 
been solved long ago using both semi-classical and 
fully quantum-mechanical approachsi 6 i 7 ' 8 Indeed, the 
fully quantum-mechanical, relativistic case of scatter- 
ing two plane-wave electrons has long been a textbook 
problem This classic problem was revived recently in 
the works of Jouffrey et alX and of Schattschneider et 
alr^ in which a "flaw" in the standard theory is pointed 
out. The flaw is the approximation that the so-called 
"longitudinal" and "transverse" matrix elements for the 
scattering process may be summed incoherently, as ar- 
gued by Fano in a seminal paper^ In fact, this approxi- 
mation is only valid when the sample under consideration 
posseses certain symmetries. In a later review article, 8 
Fano states this condition explicitly; namely that his 
original formula for the cross-section is only applicable to 
systems of cubic symmetry. However, this caveat, seems 



to have been generally ignored, and hence turns out to be 
the source of the magic angle "mystery" X Jouffrey et al., 
and later Schattschneider et al., showed that if one cor- 
rectly sums and squares the transition matrix elements 
then, in the dipole approximation, one finds the magic 
angle corrected by a factor of two. 

Our aim here is to examine the theory in more detail in 
order to derive both relativistic and material-dependent 
corrections to the magic angle. In Section II we consider 
relativistic electron scattering within the formalism of 
quantum electrodynamics (QED). Working the Coulomb 
gauge, we show that one can almost reproduce the results 
of Jouffrey et al. and the theory of Schattschneider et al. , 
apart from a simple correction term of order hu/mc 2 , 
which is not always negligible. Here Tiw is the energy lost 
by the probe and mc 2 is the rest energy of an electron. 
In Section III we suggest the possibility of incorporating 
macroscopic electrodynamic effects into the theory, which 
can break the symmetry of sample independence of the 
magic angle. 



II. COULOMB GAUGE CALCULATION 

An appealing aspect of the formalism of Schattschnei- 
der et al. is its simplicity. Their approach is similar 
to the semi-classical approach of M0ller, 5 but with the 
added simplification of working with a probe and sam- 
ple described by the Schrddinger equation, rather than 
the Dirac equation. They also find that the theory is 
simplified by choosing to work in the Lorentz gauge. Un- 
fortunately, however, the theory of M0ller is somewhat ad 
hoc in that a classical calculation in the Lorentz gauge is 
modified by replacing the product of two classical charge 
densities by the product of four different wavefunctions 
in order to obtain the transition matrix element. For 
the M0ller case this procedure is justified a posteriori 
by the fact that it reproduces the correct result, but is 
only rigorously justified by appealing to the method of 
second quantization.- M0ller's proceedure is physically 
reasonable a priori, because M0ller was interested in the 
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scattering of electrons in vacuum. However, the theory of 
Schattschneider et al., which largely mimics M0ller's the- 
ory, is less physically reasonable a priori, since the elec- 
trons are not scattering in vacuum, but are inside a solid 
which can screen the electrons. Nevertheless, since the 
discrepancy is small, the Schattschneider et al., theory is 
justified a posteriori to a lesser extent by experiment. 2 
We thus refer to the theory of Schattschneider et al. as a 
"vacuum-relativisitic theory." Consequently, in an effort 
to account for the discrepancy with experiment, we feel 
that it is useful to rederive the results of Jouffrey et al. 
from a more fundamental starting point. 

It is easy to see that the theory of Schattschneider et 
al. is not formally exact, though for many materials the 
error in the vacuum relativistic limit is negligible. In 
fact, the discrepancy can be easily explained via single- 
particle quantum mechanics: although Schattschneider 
et al. work explicitly in the Lorentz gauge, they also make 
the assumption that the momentum and the vector po- 
tential commute, 



p ■ A(r) = A(r) ■ p. 



(1) 



Of course, this commutation relation is only exact in the 
Coulomb gauge. In the end, however, the error in this 
approximation only effects the final results (e.g., matrix 
elements) by a correction of order hu/mc 2 compared to 
unity, where Tilo is the energy lost by the probe. Since 
%uj I mc 2 is at most (Za) 2 for deep-core energy loss, the ef- 
fect is usually negligible, except of course, for very heavy 
atoms. To see how corrections such as the above en- 
ter into the theory, and further to determine whether 
or not such corrections are meaningful or simply arti- 
facts of the various approximations used in the theory of 
Schattschneider et al, we find it useful to present a fully 
quantum-mechanical, relativistic many-body treatment 
along the lines of FanojS but without any assumption 
of symmetry of the sample. Our treatment is at least 
as general as that of Schattschneider et al. as far as the 
symmetry of the sample is concerned. Thus going beyond 
the formulations of Schattschneider et al. and M0ller, we 
take as our starting point the many-particle QED Hamil- 
tonian. We then show that in a single-particle approxi- 
mation the theory yields the result of Schattschneider et 
al. together with the correction mentioned above. 

Our starting point therefore is the Hamiltonian in 
Coulomb gauged 



H = H C 



H T 



ad 



(2) 



where the Hamiltonian has been split into three parts: i) 
the unperturbed electron part 



H c 



d 3 x^(x) (cat ■ p + /3mc 2 ) ip(x), (3) 



where ip(x) is the second-quantized Dirac field, oti and 
(3 are the usual Dirac matrices, m is the electron mass, 
and c is the speed of light; 



ii) the unperturbed (transverse) radiation part 



k t=l 



(4) 



where ak.i destroys a photon of momentum k, polariza- 
tion £k,i, and energy hoik', and 
iii) the interaction part 



Hi nt — 



• J d 3 x ij^(x)a. ■ A(x)ip(x) 



- I d 3 xd 3 y 



\x - y\ 



(5) 



where 



A(x) = 



E hnhc 2 ( 



(6) 

e = |e| is the charge of the proton, and V is the system 
volume. 

Let us next specialize to the case of a fixed number 
(N + 1) of electrons where the (N + l)-th electron is sin- 
gled out as the "fast probe" traveling with velocity v , 
and the remaining N electrons make up the sample. We 
also introduce a lattice or cluster of ion-cores (below we 
consider only elemental solids of atomic number Z but 
the generalization to more complex systems is obvious) 
which is treated classically, and which gives rise to a po- 
tential v c . COTC (x) — "Y^J-i {— Z e 2 ) / \x — Rj\ as seen by the 
electrons. In this case our Hamiltonian becomes: 



H = \ca ■ (p + -A(r)) + (3mi 

N 

+ J2 [ ca(i) ■ (P W + - c Mr (l) )) + P 



(r) 



(''mc 2 



-0)1 



e ^2 \ r _ r (i)| + 2 51 r"! - /• 
i=l 1 1 i=i/j=i 1 

N 

^ ^c-core(^ ) ~t~ ^core-core -^rad 3 (7) 



where the coordinates which are not labelled by an index 
refer to the probe electron. The interaction ^core-core be- 
tween ion cores is a constant and is henceforth dropped. 

To proceed to a single-particle approximation for the 
sample, the interaction of the sample electrons among 
themselves and with the potential of the ion cores may be 
taken into account by introducing a single-particle self- 
consistent potential v(x) which includes both v c . COTe (x) 
and exchange-correlation effects. The interaction of the 
probe electron with the effective single electron of the 
sample will be considered explicitly. The difference be- 
tween this interaction and the actual interaction between 
the probe and sample can be accounted for by introduc- 
ing another potential v'{x) which is not necessarily the 



3 



same as v(x); v'(x) is, in theory, "closer" to the pure 
v c -corc( x ) potential than v(x) though, in practice, this 
difference may not be of interest (see the Appendix for 
further explanation of this point). The potential v'(x) 
leads to diffraction of the probe electron, which will not 
be considered in this paper in order to make contact with 
the theory of Schattschneider et al. It is also for this rea- 
son that we have introduced a single-particle picture of 
the sample, along with the fact that we want to apply 
this theory to real condensed matter systems in a practi- 
cal way. The extension to the many-body case, in which 
the only single-body potential seen by the probe is due 
to the ion-cores, is given in the Appendix. Thus using 
the single-particle approximation for the sample, 



H = 



col ■ [p + -A(r)] + /3mc 2 
c 



ca s ■ [p s + -A(r s )} + (3 s mc 2 
c 

v'(r) + e 2 - — - — , + v(r s ) + H xad , 
\r - r s \ 



(8) 



where the quantities labeled by the letter s refer to the 
sample electron and the unlabeled quantites refer to the 
probe electron. In the remainder of this paper we set 
v' — ► 0, though the generalization of the theory to include 
diffraction is not expected to be difficult. 

As it turns outji 2 - we may start from an effective 
Schrodinger treatment of both the sample and the probe 
rather than a Dirac treatment. The treatment of the 
probe by a "relativistically corrected" Schrodinger equa- 
tion is standard practice^ in much of EELS theory, and 
is appropriate^ for modern microscope energies of in- 
terest here (e.g., a few hundred keV). The relativistic 
correction to the Schrodinger equation of the probe con- 
sists in simply replacing the mass of the pr obe m by 
the relativistic mass ml = m/y where 7 = 1/yl — Vq/c 2 . 
Moreover working with a Schrodinger equation treatment 
facilitates contact with the "vacuum-relativistic" magic- 
angle theory of Schattschneider et al. We will indicate 
later how the results change if we retain a full Dirac 
treatment of the electrons. Thus we may start with the 
Hamiltonian 



H 



{p+(e/c)A(r)} 2 , [p s + {e/c)A{r s )f 



2m' 



2m 



+ v{r s ) + 



a< 1 



#0 + ~^-P ■ A{r) + — p s ■ A(r s ) 

m c mc 

, e " , + 0(A 2 ) . 



(9) 



In this theory the unperturbed states are then direct 
products of unperturbed sample electron states (which in 
calculations can be described, for example, by the com- 
puter code FEFF8,W) unperturbed probe electron states 
(plane- waves, ignoring diffraction), and the free (trans- 
verse) photon states. Also, from now on we ignore the 



interaction terms which are 0(A 2 ). Thus our perturba- 
tion is 



U = 



^— p ■ A(r) + — p s ■ A(r s ), 
> > ' 1 • 'lit /■ 



and we are interested in matrix elements of 

U + UG U+... 
where the one-particle Green's function is 

G (E) 



1 



E- H + ir) 



(10) 



(11) 



(12) 



and r\ is a postive infinitesimal. The matrix elements are 
taken between initial and final states (ordered as: probe, 
sample, photon) 



|J) = \i) |0) and \F) = \k F ) \ f) |0) 



(13) 



To lowest order (e 2 ) there will be a "longitudinal" (in- 
stantaneous Coulomb) contribution to the matrix ele- 
ment, and a "transverse" (photon mediated) contribu- 
tion, as illustrated in Fig. 1. 

Instead of elaborating the details from standard per- 
turbation theory, we simply write down the result for the 
matrix element 



M 



4ne 2 



1 



</|e^K> 



y 

n 1 



or 



2 2 

c z q z m 



7 (/I -e* 9 



m 



(14) 



where fc-r (see Fig. 2) is the part of the initial (or fi- 
nal) momentum which is perpendicular to the momen- 
tum transfer hq. In the remainder of this paper we will 
choose our units such that h = 1. 



~k? T = 



Q 2 J 



q 2 



(15) 



The result of Eq. IT4]) is easy to understand diagramat- 
ically. For example, to each wiggly line of momentum q 
and energy w we may assign a value 



n2 



2nc 



(16) 



At this point we note that the relativistic many-body 
version of Eq. Q14p can be obtained by making intuitively 
reasonable replacements such as p/m — > ca, e lq rs — > 

J2i e l<J " r<,) . See the Appendix for further details. 

Eq. (|14p is equivalent to the matrix elements given by 
Fano in Eq. (12) of Ref. 8[. The cross-section given by 
Fano in Eq. (16) of Ref. [8|, in which the matrix elements 
have been summed incoherently, is not generally correct 
and is the source of the magic angle "mystery" 




FIG. 2: The relevant momenta: ki is the initial momentum 
of the probe particle, k,F is the final momentum of the probe 
particle, q is the momentum transfer fcj — kp and fey is the 
part of both the initial of final momenta which is perpendic- 
ular to the momentum transfer. 




FIG. 1: Feynman diagrams for the scattering process due to 
both the instantaneous Coulomb interaction (upper) and the 
transverse photon interaction (middle, lower). The solid lines 
labeled by momenta kj and ki? represent the probe particle; 
thick solid lines labeled by the letters i and / represent the 
sample particle; the dashed line is the instantaneous Coulomb 
interaction; and the wiggly lines are transverse photons. Time 
flows to the right. 



Before continuing to the dipole approximation it is use- 
ful to rewrite Eq. (fT4|) using the definition 



j , qki 

k T = k I q — — 



(17) 



to eliminate fey in favor of ki (or equivalently Vq 



k]/m'). Making this replacement we obtain 



M 



■lire 2 



</|e 



iq r I 



q-v Q (f\q-pe iq " 



iq-r 



(f\vp ■ (p/m)e 
to 2 — c 2 q 2 



which can be rewritten as: 
47re 2 1 

M ■■ 



(18) 



V q 2 -{uJ 2 /c 2 ) 

,2 



1 



v P 



(f\e lq - r 

q 



q 2 c 2 



1 



P 



(19) 



where we have made use of q ■ Vq — uj in order to cancel 
certain terms which appear after commuting the expo- 
nential through to the far left. Also, we have removed the 
label s from the position and momentum of the sample 
electron. This change in notation will be used through- 
out the remainder of this paper. 

Eq. (|19p is the same as Eq. (6) of Schattschneider et 
al., except for an "extra" term 



(f\e lq - r (l 



q p 



(20) 



Fortunately, this term may be simplified by considering 
the commutator 



[e i «- r ,H ] = [e i - 



2m J 



pq 



g 

2m 



(21) 



where the first equals sign follows from the fact that e lq r 
commutes with everything in Hq except for the kinetic 
term of the sample electron (by its definition Hq explic- 
itly contains only local potentials). Then, using the fact 
that for any operator O, 

</| [O, H ] \i) - (/| O \i) (Ei-Ef) = (/| O \i) (-co) (22) 

we have 

■2 



(/| e^ r |t) (-co) = (/| e^(^ + £-) \i) (23) 



2m' 



5 



and thus 



A. Dipole Approximation and the Magic Angle 



(/ | e ^ (1 _^ )N) = (/ | e ^ ? {i) . 

rauj Zrabj 

Making the above replacement in Eq. (fTi?]) we find 
47re 2 1 



(24) 



M 



V q 2 — uj 2 / c 2 
x (/| e Wl_ 



me 2 q 2 c 2 2mu) J 



i) (25) 



and we see that the "extra" term only changes the re- 
sult by order uj/mc 2 where mc 2 is the rest energy of an 
electron and u is the energy lost; 



M 



Aire 2 



V q 2 — lu 2 /c 2 



(f\e' qr (l 



Vq p 



mc 2 2mc 2 ) 



Aire 2 



V q 2 -u 2 /6 



Kr^{f\e iqr (l-^-(p+h)\i). (26) 
ur/c \ mcr 2 / 



Eq. ([26)1 is the same as what Schattschneider et al. would 
have obtained if they had not neglected the commutator 
[p,A]. 

That a term proportional to p + q/2 rather than sim- 
ply p appears in Eq. (|26|) is correct and can be under- 
stood from the following simple example: The interaction 
Hamiltonian for a point particle with an external field is 
given by e<f> — ev ■ A/c, or rather 

#int ~ J d 3 x (n(x)ct>(x) - - c j{x) ■ A{x)^j , (27) 

where n(x) is the density and j{x) is the current, and 
where the above integral, with the potentials considered 
as functions of the source location, is a convolution in 
space and thus a product in Fourier space-the rough cor- 
respondence indicated by the "~" symbol in Eq. (|27p is 
considered more rigourously in the Appendix. Next, we 
note that the Fourier transform of the current density (in 
second-quantization) is given for a free particle by^ 



i 

mV 



Cfe 



where 



4>(x) = y^Cfct 



ikx 



(28) 



(29) 



and where q is considered to be the momentum trans- 
ferred to the sample. This is in agreement with the usual 
conventions of EELS 



kr — ki 



(30) 



Thus we see that Eq. (|26l) is indeed correct, in both sign 
and magnitude of the "extra" term. 



In the dipole approximation Eq. (|26|) reduces to 



1 



V q 2 — uj 2 /c'' 



(/I (i 



q ■ r 



P 

mc 2 



(31) 



The term (1 — Vq ■ q/2mc 2 ) does not contribute because 
(i\f) = 0. Now, we make use of the replacement p/m — > 
iojr which is appropriate within the matrix element to 
find 



47re 2 



V q 2 - uj 2 /c' 



(/I Q 



Vo(q ■ v Q ) 



r i 



(32) 



We have thus found the same "shortened g-vector" that 
appears in Eq. (15) of Schattschneider et al. and Eq. (2) 
of Jouffrey et al. Specifically, for an initial electron veloc- 
ity Vq in the z-direction, we have found the replacement 
q z — * q z (l — Vq/c 2 ) which in turn leads to a significant 
correction (on the order of 100 percent for typical elec- 
tron microscopes) to the magic angle. 

The magic angle 9m is defined for materials with a "o 
axis" by the equality of two functions of collection angle 



F(a c ) = / dee 



[e 2 + e 2 E /^] 



(33) 



and 



G(a c ) = 2 



2 ra c | 

- n " E 1 dee 



. . - . (34) 



T Jo [o 2 + e 2 E / 1 4 



where 7 = (1 — u 2 /c 2 ) -1 / 2 , 0e is the so-called "character- 
istic angle" given in terms of the energy-loss w, the initial 
probe speed vq, and k^E — u>/vq. Both of the above in- 
tegrals may easily be evaluated in terms of elementary 
functions, but we leave them in the above-form for com- 
parison with the theory of the Section III. Eqs. ([33| and 
(fM)) both make use of the approximation sin(0) e. 
Since typical scattering angles are on the order of milli- 
radians this small angle approximation is highly accurate. 

The expressions for F(a c ) and G(a c ) are easily de- 
rived within the framework of the Schattschneider "vac- 
uum theory"— and result in a ratio of magic angle to 
characteristic angle which is independent of the material 
which makes up the sample. The factors of (1 — Vq/c 2 ) 
which appear in Eqs. (|33[) and (|34[) come from including 
the transverse effects (as in Section I) and thus the non- 
relativistic (c — * 00) result for the ratio of magic angle to 
characteristic angle is independent of transverse effects. 
The "transverse" correction to the magic angle is on the 
order of 100 percent. This corrected theoretical magic 
angle is in much better agreement with the experimen- 
tally observed magic angle, although the experimentally 
observed magic angle seems be somewhat larger (on the 
order of 30 percent) and sample dependent.— These fur- 
ther discrepancies between theory and experiment are 
addressed in Section III. 



6 



III. MACROSCOPIC ELECTRODYNAMIC 
EFFECTS 

As discussed above, the result of Schattschncider et 
al. is nearly in agreement with that obtained in Section 
II of this paper in the vacuum relativistic limit. However, 
because of the residual discrepancy between these results 
and experiment we now consider how macroscopic elec- 
trodynamic effects can be incorported into the quantum 
mechanical single-particle formalism. We find that the 
corrections to the magic angle which result can be quite 
substantial at low energy-loss. However, we are unaware 
of any experimental data in this regime with which to 
compare the theory. Nevertheless, the inclusion of di- 
electric response introduces a sample dependence of the 
theoretical magic angle which is consistent with the sign 
of the observed discrepancy. 

Certain condensed matter effects are already present 
in the existing formalism via the behavior of the ini- 
tial and final single-particle states in the sample, and 
in many-electron effects which are neglected in the inde- 
pendent electron theory. However, the macroscopic re- 
sponse of the sample can be taken into account straight- 
forwardly within a dielectric formalism. This procedure 
is similar to the well-known "matching" procedure be- 
tween atomic calculations and macroscopic-dielectric cal- 
culations of the stopping powerj 16 i 17 i 18 That is, the fast 
probe may interact with many atoms at once, as long 
the condition vq » u^a (where loq is a typical elec- 
tronic frequency and a a typical length scale) is fulfilled. 
Under these conditions the sample can be treated using 
the electrodynamics of continuous medial 

Effects due to the macroscopic response of the system 
can be included within a formalism that parallels that of 
Schattschneider et al. simply by choosing the "general- 
ized Lorentz gauge"— for a given dielectric function e(w), 
instead of the Lorentz gauge of the vacuum-relativistic 
theory. In the generalized Lorentz gauge, most of the for- 
mal manipulations of Schattschneider et al. carry through 
in the same way, except that instead of Eq. (|19p we end 
up with 



M 



Aire 2 



1 



e(u))V q 2 — e(aj)to 2 /c 2 



</|e" 



qr 



(35) 



The factors of e in Eq. ([35)) can be understood physically 
as due to the fact that c — > c/y/e in the medium, and also 
to the fact that the sample responds to the electric field 
E rather than the electric displacement D. Eq. (|3"5|) is 
derived in the following subsection. 



A. Generalized Lorentz Gauge calculation 

We consider a probe electron which passes through 
a continuous medium characterized by a macroscopic 



frequency-dependent dielectric constant e(w) and mag- 
netic permeability p = 1. It is appropriate to ignore the 
spatial dispersion of the dielectric constant at this level 
of approximation. 1 - Then Maxwell's equations are 



with D — eE. And 

V x B 



V D = i%p ext 
47rj ext 



1 dD 

c~d7' 



(36) 



(37) 



where the charge/current densities p cxt and j ext refer 
only to the "external" charge and current for a probe 
electron shooting through the material at velocity Vq. 
The other two Maxwell equations refer only to E and B, 
and can be satisfied exactly using the definitions 



-V<j>- 



1 dA 
c~dt ' 



and 



B = V x A 



We next insert Eqs. fl35J| and into Eqs 
and choose the generalized Lorentz gauged 

V- A+-^- [ dt'e(t-t'U(t')=0 
cot J 



(38) 

(39) 
and (J37D 

(40) 



This gauge choice leads to the momentum space (q,u>) 
equations 



2 , / \' 



An 



and 



A(q, uj) = An 



Pcxt(g ; ^) 

Jext(g,w) 



(41) 



(42) 



We now write p ext ,(q,u) = (-2ne)S(ui- q-v ) and j ext = 
voPext to find explicit expressions for <f> and A: 

, si, \ An{—2ne)5(uj — q ■ Vq) , , 



[e(u>)uj 2 / c 2 ] — q 2 



and 



Vo 



A(Q^) = — eM</>(g,w) 



(44) 



Then, proceeding roughly in analogy with Schattschnei- 
der et al, we have 



H 



Ho 



Ha 



2mc 

e 

2mc 



(p ■ A + A ■ p) — e<p+ 0(A 2 ) 

(2A ■ p — iV • A) — e<f>+ 0{A 2 ). (45) 



Next, evaluating the perturbation U = H — H with A = 
(vo/c)e(u)4>, we find 



U = — 

mc 



e(w) . Vq 

) Vq p- It— ■ V0 

c 2c 



e<f) 



(46) 
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In calculating the matrix element of U it is appropriate 
to replace V<p by f° r the case when the final states 
are on the left in the matrix element. Thus 



M = (f\(kf\U\h)\i)=-e<t>(q,u;) 



(/!■ 



Alternatively, since 



7 v ■ (p + q/2) 



e(u>) (q 2 — u) 2 e(uj)/c) 



(47) 



(48) 



we have 



M 



Aire 2 



e(u>)(q 2 — c(uj)lu 2 / c 2 ) 



(f\e 



iqr 



T V ■ (p + q/2) 



(49) 



In the dipole approximation Eq. (|49|) reduces to 



4ne 2 



1 



e{u)V q 2 - e{oj)uj 2 /c' 



(f\i 



e(u) 



vq(q ■ «o) 



(50) 

where e(uj) is the generally complex valued macroscopic 
dielectric constant as which can be calculated, for exam- 
ple, by the FEFFOP— code. Consequently we find that 
that instead of the longitudinal g-vector replacement 



found by Jouffrey et al 
obtain the replacement 



>q z (l-(3 2 ) (51) 
and Schattschncider et al., we 



q z [l-e{u)f3 2 



(52) 



which is appropriate for an electron traversing a contin- 
uous dielectric medium. In the same way that Eq. (|5Tj) 
can be understood classically as being due to a charge in 
uniform motion in vacuum^ Eq. (|52p can be understood 
as due to a charge is in uniform motion in a medium. Be- 
cause the motion is uniform, the time dependence can be 
eliminated in favor of a spacial derivative opposite to the 
direction of motion and multiplied by the speed of the 
particle. For motion in the z-direction 



d d 
ot az 

Therefore, if we consider the electric field 
1 dA 



E = -V0- 



dt 



-V0- 



v dA 

c dz ' 



(53) 



(54) 



Eq. ([STj) follows from the substitution A z — (vo/c)(f>, 
whereas Eq. (f52"|) follows by making the correct substi- 
tution in the presence of a medium 



(55) 



which in Fourier space gives 



E{q,Lo) 



iq + z 



Auj)v 2 



iq. 



(q,Lu), (56) 



which is equivalent to Eq. ([521) . 

Because Eq. (|52|) depends on the macroscopic dielectric 
function the ratio Om/Qe, which formerly was a function 
only of vq, will now show material dependence. This 
is seen from the generalization of Eqs. (|3"3"j) and (|3"4")> . 
the equality of which gives the magic angle. Instead of 
Eq. (|3"3"]) for F(a c ) we now have 



dee- 



and, instead of Eq. (|34|) we now have 



G(a c ) = 29%\g\' 



\9 2 



where 



g=l- e(uj)vl/c 2 



(57) 



(58) 



(59) 



is a complex number which replaces I/7 2 in the vacuum 
relativistic theory. 

If one can calculate the macroscopic dielectric func- 
tion of the sample by some means^ then the material 
dependent magic angle can be determined theoretically 
and compared to experiment. Furthermore, the correc- 
tion to the magic angle given by the introduction of the 
macroscopic dielectric constant relative to the relativis- 
tic macroscopic "vacuum value" of Jouffrey et al. is seen 
to be typically positive (since Re[e] < 1 and < Im [e]), 
in rough agreement with observation^ In fact, it turns 
out that the correction is always positive for the materi- 
als we consider and is substantial only for low energy-loss 
where the dielectric function differs substantially from its 
vacuum value. For modern EELS experiments which use 
relativistic microscope energies and examine low energy- 
loss regions, the effect of the dielectric correction on the 
magic angle should be large. 

Example calculations using our relativistic dielectric 
theory compared to both the relativistic vacuum theory 
of Schattschncider et al. and to the non-relativistic vac- 
uum theory are shown in Fig. ^ for the materials boron 
nitride and graphite. The data of Daniels et al~ is also 
shown in the figures. We have not attempted to estimate 
the true error bars for the data; the error bars in the fig- 
ure indicate only the error resulting from the unspecified 
finite convergence angle. 



IV. CONCLUSIONS 

We have developed a fully relativistic theory of the 
magic angle in electron energy loss spectra starting from 
the QED Hamiltonian of the many body system. As 
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with the single-particle theory of Jouffrey et al. and 
Schattschneider et al. we find a factor of two "transverse" 
correction to the non-relativistic ratio 0m/&e- We have 
also shown how macroscopic electrodynamic effects can 
be incorporated into the relativistic single-particle for- 
malism of Schattschneider et al. In particular we predict 
that these dielectric effects can be important for deter- 
mining the correct material-dependent magic angle at low 
energy-loss, where the difference between the dielectric 
function relative to its vacuum value is observed to be 
substantial. 

Several other factors may be important for correctly 
describing the energy loss dependence of the magic an- 
gle in anisotropic materials. In particular, we believe 
that further study of the many body effects (beyond a 
simple macroscopic dielectric model) via explicit calcula- 
tions of the microscopic dielectric function and including 
time-dependant density functional/Bethe-Salpeter the- 
ory TDDFT/BSE 22 are an important next step in the 
description of all EELS phenomena, including the magic 
angle. 



2 
1.5 
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50 100 150 200 250 

Energy-Loss (eV) 

FIG. 3: The magic angle to characteristic angle ratio 0m/6e 
is compared for three differing theories and one experiment. 2 
The materials considered in the figure are boron nitride (top 
figure) and graphite (bottom). The microscope voltage is 
fixed at 195 keV. Both the non-relativistic and relativistic 
vacuum theories show no dependence on the energy-loss and 
no dependence on the material. The relativistic dielectric 
theory shows that the magic angle should deviate from the 
vacuum value by a significant amount in regions where the 
macroscopic dielectric response is substantial. 
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APPENDIX: RELATIVISTIC EFFECTS 

Starting from Eq. ([7]) we write (the notation Ho in the 
Appendix differs from that in the main body of the text): 

2 i 



H 



ca.-p + f3mc + v c . COIC (r) 

N 



- E 

2 V 



and 



| r (i) _ r (i)| 



N 



ad • 



(A.l) 



N 



U = ea- A(r) + ^ 



ea (>) ' A{r') 



i=i 



We are interested in matrix elements of the perturbation 

U + UG Q U + ... (A.3) 

between eigenstates of the unperturbed Hamiltonian 

\I) = \h) |*<) |0) and \F) = \k f ) \9 f ) |0) . (A.4) 

The difference between the many-body case and the 
single-particle theory of the sample is that the wavefunc- 
tion of the sample now depends on N electron coordi- 
nates, instead of one effective coordinate. Also we see 
that the only potential "seen" by the probe (i.e., included 
in the unperturbed probe Hamiltonian) is the w e - core po- 
tential. This is to be contrasted with the "unperturbed" 
sample Hamiltonian which includes not only the fc-corc 
but also the Coulomb interactions between all the sample 
electrons. 

Consequently, working with a unit volume and pro- 
ceeding exactly as in the single-particle case, we find a 
"longitudinal" contribution to the matrix clement 



Mr 



Aire 2 



N 



i=i 



and a "transverse" contribution 
Aire 2 



N 



M T = 



uj 2 /c 2 — q 2 



(A.6) 



i=l 



where 



olt = at — q- 



q ■ a 



(A.7) 



and where the u(p) are the usual free-particle Dirac 
spinors, normalized such that 



u'(p)u(p) = 1 



(A. 



The two matrix elements Ml and Mt are to be 
summed and then squared, but before proceeding with 
this plan we make the following useful definitions: 
The transverse Kronecker delta function (transverse to 
momentum-transfer) 



8$ = 8 lJ 



q l q> 



the (Fourier transformed) density operator 



A' 



»(«) = £« 



and the (Fourier transformed) current operator 



3(Q) = 



N 



(A.9) 



(A.10) 



(A.H) 



Next, we recall some properies of the Dirac spinors 
u(p) and of Dirac matrices which we will presently 
find useful: i) There arc four independent spinors 
{u^\ u^ 2 \ , u*- 4 -* }, the first two of which will refer to 
positive energy solutions, and the second two of which 
will refer to negative energy solutions (and are not used 
in this calculation); 

ii) the positive energy spinors satisfy a "spin sum" 



1 



s=l 



2E(p) 
1 

2E{p) 



(E(p) + col ■ p + f3mc 2 ) 
(E(p) + h D (p)) , (A.12) 



where E(p) = \Jp 2 c 2 + m 2 c 4 ; 

iii) the Dirac matrixes satisfy the trace identities 



TT(aiOCj) 



4<5 4 



(A.13) 



Tr (aiCtjaikan) = 4 (SijSki - S ik 5ji + Su5 jk ) ,(A.14) 
TT{oi r a. j ) = Tr(a^«4) = 45*?, (A.15) 
Tr(a ,; a J T a fc a^) = 4 U*8$ ~ 5^5^ + ^^ fe )(A.16) 

iv) Finally, we note that in this calculation there are 
many simplifications due to the fact that lu << mc 2 < 
E(ki) ps E(kp). For example, 



1 



2E i E j 



{EiEj+m 2 c 4 + c 2 ki ■ kj) 



E(h 



0( 



r)«i; 
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throughout the calculation we ignore terms of order 
u>/E(pi) Using these identities it is easy to see that 



\ £ £ l^l 2 = 



Si =l Sf = l 



Aire* 



(*F|nt|^)| (A.17) 



which has the same form as in the non-relativistic case 
(up to order lo/ E{)\ the squared matrix element is much 
simplified by the sum over final probe-spin and average 
over initial probe-spin. Of course, the matrix element 
itself is completely general in terms of probe-spin, but 
many simplification arise from ignoring the probe-spin 
and exploiting the spin-sums. 

Continuing on to the transverse matrix element-and 
including a few more of the details (a, b, c, and d are 
Dirac indices)-we find 



22 1 2 2 

j££i^i 2 = j££ 



s 4 =l s f = l 



Aire 2 



i s f) „,mab 



xu(k f y a 



Aire 2 



k lu 2 /c 2 — q 2 

xTr ((E(kf) + h D (k f ))a^(E(ki) + h D (k t ))a^) 
Aire 2 



^u) 2 /c 2 — q 2 
For the cross term we find 

2 2 

2 



(A.18) 



2 2 

\ £ JlfcJU? = ( 47 re 2 



1 



Si = l S/ = l 



q 2 (u 2 /c 2 — q 2 ) 



(A.19) 



The last equality follows from 

q- (Vl\j(q)\*F) = u(*i\n(q) \$ F ) , 



(A.21) 



which itself follows by considering the commutator anal- 
ogous to that of Eq. (f2"l"T). 

The final line of Eq. (|A.20[) is quite pleasing since we 
have found that if we can "ignore" the spin of the probe 
particle, we may as well have started by taking matrix 
elements between electronic states only of the much sim- 
pler interaction Hamiltonian 



U' = / d 3 x 



n(x)(j} UJ (x - x p ) 



j(q) ■ A u (x - x p ) 



(A.22) 

where the fields {<f>ui, A^} are just the e~ luJt components 
of the classical field of a point charge of velocity Vq in 
the Lorentz gauge, and where 



n(x) 



N 

£ 



S(x 



.(01 



and 



j( x ) =J2ca^5(x-x^) 



(A.23) 



(A.24) 



That is, if we take Eq. (|A.22[) as our starting point and 
proceed in the usual way, we will find that our squared 
matrix elements are exactly the same as what we know to 
be correct from Eq. (j A.20|) . The photons have dropped 
out entirely! 



Thus we have finally derived an expression for the rel- 
ativistic many-body summed-then-squared matrix ele- 
ments summed and averaged over spins, 



|££|m l + m t | 



/ 47re 2 



V v 



n(q) \Vf)_ (*/| v t ■ j(q) \V F ) 



ine 2 



lu 2 — q 2 c 2 



V{w 2 /c 2 ^q 2 ) 



(A.20) 



